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ABSTRACT 

We present striking results from Rossi X-ray Timing Explorer (RXTE) observations 
of the 21-hr low mass X-ray binary X1624— 490, showing five complex dips in unprece- 
dented detail. For the first time, dipping is detected up to 15 keV. Prominent flares are 
also observed in the light curves, limited to energies above ~8 keV. Spectra selected by 
intensity during dip episodes can be well fit with a two-component model consisting of 
a point-like blackbody from the neutron star and progressive covering of an extended 
Comptonized region, presumably an accretion disk corona (ADC), corrected for photons 
scattered into and out of the X-ray beam by a interstellar dust halo. We find that the 
outer regions of the absorber are highly ionized and that electron scattering is totally re- 
sponsible for the X-ray attenuation during shallow dipping. The timescales of dip ingress 
and egress indicate that the envelope of material absorbing the ADC has smaller angular 
size than the ADC itself, and that the ADC is likely limited to a height-to-radius ratio 
of 10%, rather than being spherical in extent. In addition, we have analyzed ~4.5 yrs 
of RXTE All Sky Monitor (ASM) coverage to derive the first accurate orbital ephemeris 
for X1624— 490, with phase zero (the time of dip centers) well-described by the relation 
245 0088.63918(69) + N x 0.869907(12) (JD). 

Subject headings: accretion, accretion disks — scattering — stars: binaries: close - 
stars: circumstellar matter — stars: individual (X1624— 490) — stars: neutron — inter- 
stellar medium: extinction — X-rays: stars 
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1. Introduction 

The study of X-ray dipping sources has played a 
critical role in the understanding and visualization of 
how a low mass X-ray binary (LMXB) might look 
to a nearby astronautical observer; it is generally ac- 
cepted that the dips are due to occultations of the 
central source by a thickened region of the accretion 
disk rim where the gas stream from the companion 
impacts upon the outer disk. However, the dips them- 
selves also provide a wealth of physical information 
not limited to their role in defining the system ge- 
ometry. From their irregular, ragged appearance and 
cycle-to-cycle differences in width and depth we know 
that the disk structure is dynamic. Also, spectral fit- 
ting results show that dips are associated with, but 
cannot be completely described by, increases in ab- 
sorption of N H k, 10 23 cm~ 2 . 

X1624— 490 is one of the most unusual members 
of this class; its persistent emission is the brightest 
(6xl0 37 erg s _1 ), its dip profiles seemingly the most 
erratic, and its 21-hr orbital period is the longest - 
~5-25 times longer than other dipping sources, cor- 
responding to a much greater stellar separation and a 
larger accretion disk radius. Dipping is deep, ~ 75% 
in the 1-10 keV band, and can at times reach a stable 
lower level (Watson et al. 1985; Church & Balucihska- 
Church 1995) suggesting that one spectral component 
is entirely removed. The source also exhibits strong 
flaring in which the X-ray flux can increase by 30% 
or more over timescales of a few thousand seconds. 
In this regard, it may provide an interesting overlap 
between the dipping LMXBs, generally assumed to 
be atoll sources by virtue of their persistent emission 
level and noise characteristics, and the Z-track sources 
with their well-known flaring branches (e.g. Hasinger 
& van dcr Klis 1989). 

The depth, duration, and spectral evolution in dip- 
ping vary considerably between members of the dip- 
ping class. An additional complication is the fact 
that several sources exhibit an unabsorbed compo- 
nent of the spectrum, e.g. XB1916— 053 (Smale et al. 
1988, 1992), XB0748-676 (Parmar et al. 1986) and 
XB1254-690 (Courvoisier et al. 1986). More recently, 
it has been demonstrated that a unified approach de- 
scribes the dipping sources (Church & Baluciriska- 
Church 1995; Church et al. 1997, Balucihska-Church 
et al. 1999). Their model consists of point-like black- 
body emission from the neutron star, plus extended 
Comptonized emission from the Accretion Disk Corona 



(ADC), and a "progressive covering" description of 
absorption. In this description, the absorber on the 
outer edge of the accretion disk (White & Swank 
1982) moves progressively across the emission regions, 
so that at any stage of dipping, part of the Comp- 
tonized emission is absorbed and part is not, giving 
rise to the observed unabsorbed component. As a 
point source, the blackbody is rapidly absorbed once 
the envelope of the absorber reaches the neutron star. 

In previous observations of X1624— 490 with EX- 
OSAT and Ginga, a blackbody plus bremsstrahlung 
model was used to parameterize the spectra (Jones & 
Watson 1989). Subsequently, Church & Balucihska- 
Church (1995) showed that the light curve at energies 
>5 keV was dominated by flaring, which can strongly 
modify the spectrum and make the spectral investiga- 
tion of dipping difficult. By selecting sections of non- 
dip and dip data without apparent flaring, the above 
two-component model was found to fit the data well, 
showing that in deep dipping the blackbody compo- 
nent was totally absorbed, and the Comptonized com- 
ponent relatively little absorbed. However, with only 
one non-dip and one deep dip spectrum it was not 
possible to determine the extent of absorption of the 
Comptonized component. 

A more recent BcppoSAX observation revealed the 
dust scattered halo of the source demonstrated by An- 
gelini et al. (1997). Modeling of the radial distribution 
of count rate in the MECS at a series of energies below 
5 keV allowed an optical depth to dust scattering of 
2.4 ± 0.4 at 1 keV to be derived, the data being consis- 
tent with a cross-section varying as E~ 2 (Balucinska- 
Church et al. 2000a). Spectral models of non-dip and 
dip data included the effects of this interstellar dust 
scattering, and good fits were obtained with the above 
two-component model, in which the blackbody had 
temperature kT — 1.31 ±0.07 keV, and the Comp- 
tonized emission was approximated by a cut-off power 
law with photon index of T = 2.01qy an< ^ cut-off en- 
ergy ~12 keV. While more sensitive than previous 
studies, the BcppoSAX observation also included only 
a single dip episode. 

Previous observations of X1624— 490 have been 
limited either by the short duration of the observa- 
tions or by the signal-to-noise ratio of the data; a 
thorough investigation requires a baseline of several 
orbits, coupled with a collecting area and efficiency 
large enough to study the spectral variability on short 
timescales. In this paper, we present light curves ob- 
tained with RXTE covering 4.5 contiguous orbits of 
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X1624— 490, during which both strong dipping and 
strong flaring took place. We derive an accurate sys- 
tem ephemeris, and perform a detailed study of the 
spectral variability in dips. A companion paper will 
contain our analysis of the evolution and physics of 
the flares (Balucihska-Church et al. 2000b, hereafter 
Paper II). 
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Fig. 1.— The X-ray light curve of X1624-490, 1999 
September. The energy range covered is 2-20 keV, 
and the time binning is 64 sec. 



2. Observations 

We observed X1624-490 with RXTE (Bradt, Roth- 
schild, & Swank, 1993) from 1999 September 27 18:22 
UT - October 1 13:11 UT, for an on-source total good 
time of 212 ksec. The X-ray data presented here were 
obtained using the PCA instrument with the Stan- 
dard 2 and Good Xenon configurations, with time res- 
olutions of 16 sec and < l^sec respectively. The PCA 
consists of five Xe proportional counter units (PCUs) 
numbered PCUs through 4, with a combined to- 
tal effective area of about 6500 cm 2 (Jahoda et al 
1996). Only PCUs and 2 were reliably on through- 
out our observations, and we limit our analysis to data 
from these detectors. Data were also obtained with 
the HEXTE phoswich detectors, which are sensitive 
over the energy range 15-250 keV (Rothschild et al. 
1998). However, because of the cut-off spectral shape 
of X1624— 490, the counts detected were limited to 
energies <25 keV, and the HEXTE data did not in 
this case provide any information additional to that 
gained using the PCA. 



Based on a preliminary determination of the source 
ephemeris, derived from the above observations and 
the RXTE All Sky Monitor (ASM) light curve, we 
scheduled the balance of our RXTE pointed observing 
time to capture a complete dip episode, uninterrupted 
by Earth occultations and South Atlantic Anomaly 
passages. These observations were performed on 1999 
November 12:20-16:21 UT, for a total good time of 
14.3 ksec. For this observation, PCUs 0, 2, and 3 were 
reliably on throughout. 
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Fig. 2. — The 1999 November light curve in the 
2-20 keV band, with 16 sec time resolution, show- 
ing the wealth of detail in a single dip episode from 
X1624-490. 



We performed our data analysis using the RXTE 
standard analysis software, FTOOLS 5.0. Back- 
ground subtraction of the PCA data was performed 
utilizing the "skyvle/skyactiv 20000131" models gen- 
erated by the RXTE PCA team. The quality of 
the background subtraction was checked in two ways: 
(i) by comparing the source and background spec- 
tra and light curves at high energies (50-100 keV) 
where X1624— 490 itself no longer contributes de- 
tectable events; and (ii) by using the same models to 
background-subtract the data obtained during slews 
to and from the source. We conclude that our back- 
ground subtractions in the 2-20 keV energy range are 
accurate to a fraction of a count per second. Since 
one of our goals was a sensitive spectral analysis, we 
took pains to check the quality of both the back- 
ground subtraction and the response matrices used, 
including a careful analysis of Crab calibration obser- 
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vations on 1999 September 26, October 13, November 
8, and November 23, bracketing our observations of 
X1624-490. Based on this work, we added 2% sys- 
tematic errors to the spectral fits. 

Our orbital ephemeris calculations made use of 
two standard tools for establishing periodicities: a 
discrete Fourier transform technique optimized for 
unevenly-sampled time series (Scargle 1989), and a 
period-folding code utilizing the L-statistic (Davies 
1990). Errors were calculated using the formulas in- 
trinsic to these methods, and corroborated by fitting 
periodic functions to the unfolded light curve data 
and calculating the formal errors on the orbital pe- 
riod and epoch directly. 

The search for fast variability in X1624— 490 uti- 
lized the high time resolution PCA data obtained in 
the Good Xenon mode. Discrete power spectral den- 
sity distributions (PSDs) were calculated by dividing 
the data into segments of uniform length, perform- 
ing fast Fourier transforms of each, and averaging the 
results. The PSDs were normalized such that their 
integral gives the squared RMS fractional variability 
(Miyamoto et al 1991; van dcr Klis 1989). We sub- 
tracted the Poisson noise level from the power spectra, 
taking into account the modifications expected from 
PCA detector deadtime. 

3. Results 

3.1. Overall source behavior 

In Figure 1, we show the background-subtracted 
2-20 keV light curve for the 1999 September observa- 
tion. Superimposed on the mean persistent flux lev- 
els, the most prominent features of the observation 
are four regular dips separated from each other by 
^21 hrs, representing the orbital period of the system, 
and strong irregular flaring, most evident at phases 
^=0.30-0.75 (assigning (f>=0 to the dip center). The 
1999 November observation is shown in Figure 2, and 
provides a good illustration of the prodigious amount 
of detail in an individual dip from this source. The 
complex phenomenology seen in the overall light curve 
becomes more amenable to examination when we di- 
vide the 1999 September data into four energy bands 
covering 2-6, 6-10, 10-15, and 15-20 keV (Figure 3). 
(This is, incidentally, the first time the variability of 
the source above 10 keV has been accessible with good 
signal to noise.) From this Figure, we can come to the 
following qualitative conclusions: 

• Dips are deepest at lower energies, but are still 



(barely) discernable above 15 keV. 

• The flaring is a higher-energy effect, not visible 
in the 2-6 keV range, but becoming the dominant 
feature in the 10-15 keV and 15-20 keV light curves. 
Examining light curves with even finer energy reso- 
lution, we can determine that flaring first becomes 
significant above 8 keV. 
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Fig. 3.— The X1624-490 light curve from 1999 
September, broken down into several bands to show 
the energy dependence of the dipping and flaring ac- 
tivity. 



3.2. The ephemeris of X1624 490 

Applying the period search algorithms to the 1999 
September PCA observation, we find a best-fitting 
period of 20.98±0.14 hrs. This is consistent with 
previous estimates from datasets of limited dura- 
tion, but for a useful long-term ephemeris, we turn 
to the publicly-available data from the RXTE ASM. 
The ASM (Levine et al. 1996) consists of three sim- 
ilar Scanning Shadow Cameras, sensitive to X-rays 
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in the 2-12 keV energy band, that cover ~80% of 
the sky every 90 minutes, with each datapoint repre- 
senting a flux measurement from a single 90s dwell. 
X1624-490 is detected in the ASM with a weighted 
mean net count rate of 3.61±0.01 count sec -1 . (For 
comparison, the Crab Nebula produces a count rate 
of ^75 ASM count sec -1 .) For this work, we analyzed 
data spanning the 1616 days from 1996 January 5 to 
2000 June 8. The 21-hr orbital periodicity shows up 
strongly even in a rudimentary Fourier transform of 
the data at a frequency of (1.3306±0.0002) x 10 -5 Hz. 
Combining Fourier and period folding analysis, we 
derive an orbital ephemeris for X1624— 490 of 

245 0088.63918(69) + NxO. 869907(12) (JD). 

This orbital period (20.8778±0.0003 hrs) is a two- 
significant-figure improvement over previous deriva- 
tions, with a predictive value limited only by the 
phase jitter and irregularity of individual dips. It is 
sufficiently accurate to extend back to e.g. the EX- 
OSAT observations of X1624-490 obtained in 1985 
March 25-27, with no cycle-count ambiguity. 

3.3. A QPO search 

Recently, a new type of quasi-periodic oscillation 
(QPO) has been discovered in three of the dipping 
sources. The feature was first seen in X1323— 619, 
where a ~1 Hz QPO was detected with an rms ampli- 
tude of 9%, remaining roughly constant in the persis- 
tent emission, the dips, and during the bursts (Jonker, 
van der Klis & Wijnands 1999). Similar features 
have been observed in X0748— 676, with a central fre- 
quency v that varied between 0.6 and 2.4 Hz and 
rms amplitudes of 8-12%, and in X1746— 370, with 
1^=1.0—1.6 Hz and rms amplitude 7% (Homan et al. 
1999; Jonker et al. 2000). The QPO in X0748-676 are 
present in all the observations analyzed by Homan et 
al. (1999) except the pointing with the highest coun- 
trate. For X1746— 370 the picture is even more clear- 
cut; the atoll-source "island" and "banana" (Hasinger 
& van der Klis 1989) branches are well defined by the 
data, and the QPOs are absent during the latter. 

For all three sources, there is little or no depen- 
dence of the QPO strength on the photon energy. 
This is very unusual: the amplitude of all other types 
of QPO seen in Z sources, atoll sources, and black 
hole candidates (BHC) from 0.01-1200 Hz shows a 
strong dependence on photon energy. Clearly the dip- 
per QPO are a different phenomenon, related to the 
high source inclination, and perhaps linked directly 



to the modulating effect of material in the accretion 
stream or at the disk edge; the exact mechanism is 
still open to conjecture. As this QPO is not observed 
in other LMXBs at lower inclinations, variations in 
the mass accretion rate are unlikely to be directly re- 
sponsible, although the disappearance of the QPO at 
higher inferred mass accretion rates may be due to 
changes in the accretion geometry. The most likely 
candidate is some modulating effect caused by the 
accretion stream, or partial covering of an extended 
X-ray source by a near-opaque medium; this struc- 
ture may disappear, change size, or change in optical 
depth at higher accretion rates. 

We have searched the X1624— 490 data for evidence 
of similar QPO activity, and find none. The extensive 
data set at our disposal allows us to place a rather 
stringent upper limit of 1% on the rms amplitude of 
a feature between 0.5-2 Hz with a FWHM of 0.4 Hz. 

The absence of such features in X1624— 490 is per- 
haps not surprising; 1 Hz QPOs are observed in 
XB1323-62 and XB1746-37 in the low-intensity is- 
land state of these atoll sources, but arc not visible 
in the higher-intensity banana state. The derived 
hardness-intensity and color-color diagrams for the 
current observations of X1624— 490 are dominated by 
the dip and flare activity and give no direct indication 
of which of the two states the source may be in, how- 
ever the high overall source luminosity (see below) 
implies a high accretion rate, perhaps incompatible 
with the production of such QPO. 

3.4. Spectral evolution during dipping 

The energy bands in Figure 3 were chosen to em- 
phasise and isolate the dipping activity, which domi- 
nates the light curve in the lowest band, and flaring, 
which dominates in the two bands above 10 keV. From 
Figures 2 and 3, it can be seen that there is a high 
degree of fast variability in dipping on timescales of 
~32 s, and that there is a tendency for dipping to sat- 
urate at a constant low level. Moreover, the 2-6 keV 
light curve shows that dipping consists of two stages: 
deep dipping with associated rapid variability, pre- 
ceded and followed by relatively shallow shoulders to 
the dipping which last ^12,000 s. Flaring, although 
apparently random in occurrence, is generally not ob- 
served during dip episodes at the high level seen out- 
side dips. However, there appears to be a section of 
flaring data at ~ 1.4 x 10 5 s, where the flaring coin- 
cides with the weak shoulder of the second dip. 
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We chose to study the spectral evolution through 
the dip episodes in the 1999 September observation 
in two stages. First, we accumulated spectra from a 
limited subset of the data, ruthlessly eliminating any 
data sections with even a suggestion of flare contami- 
nation. This approach produced 4 spectra, one for the 
non-dip emission, two at intermediate dipping levels, 
and a saturated dip spectrum. We fit these simulta- 
neously using a wide range of models to allow elimi- 
nation of most of these, and determine the probable 
form of the best-fit model. In the second stage, we di- 
vided the bulk of the observation into intensity bands 
and fit these jointly using the model determined from 
the first stage. This approach gave us confidence in 
our data selection criteria and enabled us to perform 
a sensitive search for subtle features. 

3.4-1- Stage 1: Simultaneous fits to 4 spectra 

Because of the strong spectral changes associated 
with flares, it was essential to avoid any possibility 
of data contamination with flare activity in our first 
pass. We thus selected data from the region around 
<~ 1.7 x 10 5 s (in Figure 1), consisting of a complete 
satellite orbit of data during which no flaring is de- 
tected. A deep dip spectrum was obtained from that 
part of the third dip which reached a saturated, stable 
lower level in the low-energy light curve. All traces 
of dip ingress and egress were removed by selecting 
320 seconds of data beginning at 2.1678 x 10 5 s from 
the start of the observation. Other dip spectra were 
selected in intensity bands at 180-190 count s -1 and 
130-140 count s _1 . The total source counts contained 
in the 4 spectra selected, from non-dip to the deep- 
est dipping respectively, were ~7.4xl0 5 , 3.5xl0 4 , 
4.4xl0 4 , and 2.2xl0 4 . 

These four spectra were first fitted simultaneously 
using simple, one-component models. In all cases, 
the parameters specifying the emission process such 
as the power law photon index were chained between 
the four spectra, since these cannot change during 
dipping. Results for all model fitting are shown in 
Table 1. (All values of the hydrogen column density 
quoted in the tables are in units of 10 22 atoms cm~ 2 , 
and errors are quoted to 90% confidence.) 

Simple models were unable to fit the spectra: an 
absorbed cut-off power law gave a x 2 per degree of 
freedom (dof) of 8104/188; an absorbed bremsstrahhmg 
model gave a x 2 /dof of 8789/189 and an absorbed 
blackbody model gave a x 2 /dof of 12666/189. Next, 
the two-component model consisting of a blackbody 



plus cut-off power law was tried, first with the non-dip 
spectrum alone. A good fit was obtained with x 2 /dof 
= 20/46, and so this model was applied to simultane- 
ous fitting of all four spectra, producing an overall fit 
statistic of x 2 /dof = 231/180. The two-component 
model has the form: AG(AB.BB + PCF.CPL), where 
AG represents Galactic absorption, BB the blackbody 
component, CPL the cut-off power law, AB addi- 
tional absorption during dipping and PCF the pro- 
gressive covering factor. This is the model shown to 
be able to describe spectral evolution in dipping in 
other LMXB dipping sources (see Introduction), and 
in this model the point-source blackbody is under- 
stood to be absorbed rapidly when the envelope of 
the absorber covers the neutron star, whereas the ex- 
tended Comptonized emission is covered progressively 
as the absorber moves across the extended source. 

Strong broad, residuals at ~6.4 keV were found, at 
the level of 4% in the non-dip spectrum, too high to 
be due to residual uncertainty in the instrument re- 
sponse function, and so a Gaussian line was included 
to the above model as a source term. Because of 
a tendency for widths of broad lines to increase in 
spectral fitting by absorption of continuum emission 
into the line, it was necessary to fix the width, and 
tests showed that 0.4 keV was appropriate. Fitting 
this model to 4 spectra simultaneously required the 
normalization of the line to decrease systematically in 
dipping, from 2.59xl0 -3 photon cm~ 2 s _1 in non-dip 
to 5.17xl0~ 4 photon cm~ 2 s _1 in deep dipping. Fit- 
ting the deep dip with the latter normalization gave 
X 2 /dof = 38/45, whereas fixing the normalization at 
the non-dip value gave x 2 /dof = 130/46, demonstrat- 
ing the significance of the decrease. Thus, we included 
the line within the progressive covering term to give 
a model: AG(AB.BB + PCF(CPL + GAU)). Using 
this approach, no change in the line normalization 
was necessary, showing that the line was well mod- 
eled by assuming that it originates in the same region 
as the Comptonized emission, i.e. the ADC, and so is 
subject to the same covering factor. 

At this stage, the fits still exhibited excesses at 
low energy, resulting in a quality of fit that was still 
not acceptable. Consequently, we next added the 
same interstellar dust scattering terms to the spec- 
tral model that were used in the BeppoSAX study 
(see Balucihska-Church et al. 2000a, and above), in- 
cluding terms both for scattering out of the beam 
and into the beam, giving the final form of the fitted 
model: 
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AG c- r ( AB . BB + PCF ( CPL + GAU )) + 
AG ( 1 - c- T )( BB + CPL + GAU ). 
The factor e~ T represents scattering out of the beam, 
while the (1 — c~ T ) term is scattering into the beam. 
For non-dip emission, we make the usual assumption 
that the intensities of these components are equal 
(e.g. Martin 1970). (Non-standard spectral compo- 
nents e~ T and (1 — c~ T ) were produced for inclusion 
in spectral fitting with the XSPEC package.) With 
the addition of the dust scattering terms, the low en- 
ergy excesses in the spectra were removed and a good 
fit obtained to all 4 spectra simultaneously with a 
X 2 /dof = 127/179. 

While our fitted model may appear over-constrained 
at first sight, it is crucial to note that in our simul- 
taneous fits, all emission terms (blackbody kT and 
normalization, T, cutoff and normalization values) are 
fixed at their non-dip values, and the dust scattering 
parameters are fixed at their SAX MEC values. Only 
the column densities and covering fraction are free pa- 
rameters, in a joint fit to 4 spectra of differing shapes 
and over a fivefold range of intensities. 

In addition, to demonstrate that the inclusion of 
the line and dust scattering do not favor the two- 
component model over the other candidate models, 
the same terms were added to the one-component 
models fits in the most favorable way, i.e. with the 
line being subject to progressive covering to account 
for its variation in intensity during dipping. The re- 
sults of fitting all models are included in Table 1. It 
can be seen that adding the line and dust scattering 
to simple models results in some improvement, but 
the fits are still unacceptable. Although the over- 
all fit of the best-fit model to the 4 spectra was ac- 
ceptable, the spectra (non-dip to deepest dip) having 
individual values of x 2 /dof of 20/46, 74/46, 60/46 
and 38/46, the intermediate dip spectra had broad 
excesses at ~5 keV. This was investigated further in 
our next step, which was to analyze spectra from the 
whole observation, with much improved statistics at 
the intermediate levels. 

3.4-2. Stage 2: Spectra from the complete observa- 
tion 

Spectra were selected in intensity bands 20 count s _1 
wide, with a non-dip spectrum taken from the band 
300-320 count s _1 . Five dip spectra were accumulated 
corresponding to intensity bands: 280-300, 260-280, 
220-240, 100-120 and 60-80 count s"\ such that the 
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Fig. 4. — Spectral fits to the non-dip spec- 
trum (squares); the intermediate dip level (260-280 
count s _1 ) representing the deepest level in the shoul- 
ders of dipping (see text - triangles); and the deepest 
dip spectra (circles). 



non-dip spectrum plus the first two of these provide 
good coverage of the shoulders of dipping, and the re- 
maining spectra span the main, deeper parts of dips. 
We excluded a section of data 60 ksec in length at the 
beginning of the observation, in which the 15-20 keV 
light curve (Figure 3) shows unusual behavior. The 
six spectra thus derived have a considerably higher 
total count rate than the 4 spectra previously used, 
allowing the relatively small intensity changes in the 
shoulders of dipping to be investigated. We applied 
our best fit model first to these shoulders. 

In Figure 4, the unfolded spectrum of the non-dip 
level, the deepest spectrum in the shoulders of dip- 
ping, and the deepest dip spectrum of all are shown. 
The dip shoulder spectrum can be approximately fit- 
ted as a constant vertical shift downwards of the non- 
dip spectrum, indicating that electron scattering is 
the dominant absorption process, and that photoelec- 
tric absorption is negligible. From the fitting it is 
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also clear that the blackbody requires zero column 
density, for all spectra in the shoulders of dipping. 
Thus, the intensity reduction in the shoulders must 
correspond to the absorber overlapping the extended 
ADC but not the neutron star, and the outer regions 
of the absorber which are involved in this overlap 
must be highly ionized. Comparison of the non-dip 
and deepest-dip spectra also reveals a vertical shift 
between the spectra at 20 keV showing the effects 
of electron scattering. This effect, clear in intensity- 
selected data with good statistics, was not visible in 
Stage 1. Consequently, it becomes necessary to up- 
grade the PCF model previously used to a progressive 
covering model in which X-rays are attenuated by two 
processes: electron scattering and possible photoelec- 
tric absorption expressed in the form: 

(f e -N e <r Te -N H a PE + 

where N e is the electron column density, ut is the 
Thomson cross-section, and <jpe the photoelectric 
cross-section, and the contributions of scattering and 
absorption are allowed to be independent. 

Good fits were obtained with this model, as shown 
in Figures 5 and 6. The emission parameters used are 
shown in Table 2, and we present the parameters of 
the best fits to the complete set of intensity-selected 
spectra in Table 3. The values of the emission param- 
eters used for the complete set agree well within the 
errors with the values shown for the best-fit model of 
Table 1. 

Initially the fit to the intermediate dip spectrum at 
220-240 c s^ 1 was relatively poor with an excess seen 
in the spectrum at ^5 keV and x 2 /dof = 145/45. 
This spectrum was investigated further. A concern 
throughout this work was the possibility of dip spec- 
tra being contaminated by flaring, given the strong 
spectral changes we show take place in flaring in Pa- 
per II. Deep dip spectra will not be contaminated as 
flaring is associated with the neutron star, but we 
cannot rule out the possibility for the intermediate 
spectra may be. However, it appears that this excess 
has another cause. If the absorber is highly ionized, 
residual photoelectric absorption will be due to the 
ionized absorber, not neutral absorption. In the early 
stages of dipping, this is not a problem since atten- 
uation of the Comptonized ADC emission is caused 
by the large electron column density N^' PL whereas 
the absorption column N^ PL remains zero (Table 
3) until dipping becomes deep. In deep dipping, a 




Energy (keV) 



Fig. 5. — Simultaneous fits to the non-dip and five 
dip spectra from the complete 1999 September obser- 
vation, using a two-component model consisting of a 
point-like blackbody and progressive covering of an 
extended Comptonized region, plus a halo and line 
(see text for details). Data and residuals for each in- 
put spectrum are marked as follows for easier visual 
identification: non-dip - squares; spectra of progres- 
sively deeper dipping: circles, triangles, squares, cir- 
cles and triangles. 



combination of neutral absorber and electron scatter- 
ing provides satisfactory fits, as could be expected 
since the degree of ionization must decrease inside 
the absorber. Similarly, in shallow dipping the black- 
body is not absorbed, and in deep dipping totally 
absorbed, so the question of ionized absorber does 
not arise. The exception is the intermediate spec- 
trum with N% B = 20xl0 22 atom cur 2 (Table 3). 
Refitting this spectrum with an ionized absorber for 
the blackbody resulted in a much-improved fit with 
X 2 /dof = 53/44 (as shown in the Table), and a value 
of £=118 for the ionization parameter. Strictly speak- 
ing, even this model is not totally satisfactory since 
the 1 absorb model of XSPEC that was used is based 
on a power law source whereas X1624— 490 is domi- 
nated by blackbody emission. However, we believe we 
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5 10 20 

Energy (keV) 



Fig. 6. — The unfolded spectra for the fits shown 
in Figure 5, showing the contributions of the black- 
body (BB), halo, and line to the overall fits. The 
halo components (scattered into the beam) are HBB 
and HCPL for the blackbody and cut-off power law 
contributions; the line halo is included but cannot be 
seen on this scale. 



are now at the limit of what can usefully be learned 
through spectral fitting, and further refinements will 
produce diminishing returns. 

4. Discussion 

We have shown that spectral evolution during dip- 
ping in X1624— 490 is well-described by a combination 
of the point-source blackbody plus extended Comp- 
tonization emission model, with a "progressive cov- 
ering" treatment of the absorption process. To this 
basic model it was necessary to add a broad Gaussian 
line at ^6.4 keV, and terms to describe interstellar 
dust scattering both out of and into the line-of-sight. 
The blackbody temperature kT, and power law in- 
dex T, and cut-off energy Eco from these RXTE 
fits (Table 2) are in good agreement with the val- 
ues determined from the BcppoSAX fitting, of kT 



= 1.37 ±0.07 keV, T = 2.0±g"J and E CO -12 keV 
(Balucihska-Church ct al. 2000a). During the RXTE 
observations, the total luminosity in the 1-30 keV 
band was 1.47 x 10 38 erg s" 1 at a distance of 15 
kpc (Christian & Swank, 1997), approximately twice 
that seen during the BeppoSAX observation, and thus 
complete consistency of emission parameters should 
not necessarily be expected. The bolometric black- 
body luminosity was 2.3 x 10 37 erg s _1 , or 15% of 
the total luminosity, smaller than the 36% measured 
by BeppoSAX (Balucihska-Church et al. 2000a). Lu- 
minosity considerations are discussed further in Paper 
II. 

In dipping, the blackbody component is rapidly 
and completely removed once the neutron star is cov- 
ered. The clear occurrence of saturated dipping in 
most of the dips observed indicates that in addi- 
tion to the removal of the blackbody, absorption of 
the Comptonizcd emission also reaches a stable level. 
This must be due either to due to an absorber of 
smaller angular extent than the ADC by a fixed frac- 
tion, or to the presence of blobs within the absorber 
envelope such that ~20% of the absorber is transmit- 
ting in less dense regions in between blobs of higher 
density. The Gaussian line varied in dipping in a 
way that could be described by giving it the same 
covering factor as the ADC. There are two contribu- 
tions to the excess at low energies: the dust-scattered 
halo components which are shown separately in Fig- 
ure 6, and also the part of the extended Comptonizcd 
emission that is not covered at any level of dipping. 
Table 3 shows the covering fraction rising to a maxi- 
mum value of ~80% in the best-fit model, unlike other 
dipping sources such as XB1916— 053 (Church ct al. 
1997) where dipping reaches a depth of 100% at all 
energies below 10 keV and the covering fraction sim- 
ilarly reaches 100%, proving that the angular extent 
of all emission regions is less than that of the ab- 
sorber. In the present case, there may be two rea- 
sons why this does not happen: either the absorber 
is somewhat smaller in angular extent than the ADC 
extended emission, or more likely, the absorber in the 
outer accretion disk is not uniformly dense within its 
envelope, but blobby, as indicated by the strong vari- 
ability within dipping. We investigate this possibility 
further below. 

We can obtain further information on the ADC 
from the duration of dip ingress and egress. The sharp 
ingress to deep dipping, and the fast variability within 
dipping, take place on timescales of ~32 s, and this is 
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clearly associated with absorption of the point-source 
blackbody component. However, the long, shallow 
shoulders of the dips map a gradual process of dip 
ingress and egress which is consistent with the ex- 
tended ADC being covered by relatively less dense 
absorber before the absorber envelope reaches and ob- 
scures the neutron star. The duration of the ingress 
shoulder is ~ 12 ± 2 ks from the present observation, 
and we can also derive a duration of ~ 13 ± 2 ks 
from studying the EXOSAT light curve (Church & 
Balucihska-Church 1995). Some asymmetry can be 
seen in the shoulders, and we use the shorter ingress 
time corresponding to the leading edge of the bulge 
on the accretion disk, not the trailing edge, which 
would give a longer shoulder depending on the ab- 
sorber dimensions rather than the dimensions of the 
ADC. For an absorber of larger angular extent than 
the ADC, the ingress time At is the time taken by the 
leading edge to cross the ADC given by the velocity 
of material in the outer disk: 

2 7rr disk /P = d A Dc/At, 

where rdisk is the disk radius, P is the orbital period 
and dADC is the diameter of the ADC. In the event 
that the absorber has somewhat smaller angular ex- 
tent than the ADC, as is a possible interpretation of 
spectral fitting, the ingress time gives the diameter of 
the absorber, and the source region could be as much 
as 20% larger. Using a period of 20.87 hr from our 
cphemcris work above, and a mass of 1.4Mq, a value 
of 1.0 x 10 11 cm can be calculated for the accretion 
disk radius (Frank et al. 1987). Using our observed 
At, we derive from this an ADC radius of 5.0 x 10 10 
cm. The average ingress time from the EXOSAT and 
RXTE data gives 5.3 x 10 10 cm. If the absorber has 
smaller angular size than the ADC, then the above 
calculation provides the absorber radius, so that the 
ADC could be up to 20% larger than this, on the 
basis of the covering fraction only reaching ~80% in 
deepest dipping. 

We can use the above radius of the ADC to make 
some comparisons with simple ADC theory. The 
maximum radius of an ADC that can be supported 
in hydrostatic equilibrium r cq is given by the condi- 
tion that kT < GM x m/r eq where M x is the mass of 
the neutron star and m the mass of a proton. This 
provides a simple formula: 

1.6 x 10 11 M x , , 

r cq ~ — — — (cm) 

Iadc Mq 



where Tadc is the temperature of the ADC in units 
of 10 7 K. We can use the Comptonization cut-off en- 
ergy Eco obtained in the present analysis of 12 keV 
to place limits on the electron temperature kT e of 4- 
12 keV, corresponding to an optical depth to electron 
scattering t <1 or r >1 respectively. For low opti- 
cal depth, kT e ~ Eco, but for higher optical depth, 
even if Comptonization is not saturated, kT e will be 
a factor of 2.5-3.0 smaller than Eco (e.g. Dove et al. 
1997). Using these limits on electron temperature, 
we find r oq = 1.6 — 5.3 x 10 10 cm, the smaller value 
corresponding to the higher temperature. Comparing 
these with tabc ~ 5 x 10 10 cm from our dip ingress 
calculation, we can see that there is better agreement 
with r cq calculated assuming a higher optical depth to 
electron scattering. We thus have evidence that the 
ADC is maintained in hydrostatic equilibrium. The 
ADC is likely to be "thin", with a height-to-radius 
ratio of the order of 10% at its outer edge, as opposed 
to being spherical. For example, for an ADC radius 
of 5 x 10 10 cm, it is unlikely that the absorber on the 
outer disk would be able to cover a spherical ADC 
extending to this distance in the vertical direction. 

This paper utilizes RXTE All Sky Monitor results 
made publicly available by the ASM/ RXTE Team, in- 
cluding members at MIT and NASA/Goddard Space 
Flight Center, and also archival data obtained through 
the High Energy Astrophysics Science Archive Re- 
search Center Online Service, also at NASA/GSFC. 
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Table 1 

Results from simultaneous fits to the non-dip spectrum and 3 dip spectra. 



Model 

Blackbody 
Bremsstrahlung 
Cut-off power law 
Two-cpt (prog, cov.) 

Blackbody, line, dust 
Bremsstrahlung, line, dust 
Cut-off power law, line, dust 
Two-cpt, line, dust (prog, cov.) 



N. 



H 



+0.1 

-0.0 



0.0 

4.48±0.26 
12 

lz -°-0.3 

11.6±0.7 1.55±0.04 



0+ 01 
6.0±0.1 
12.5±0.1 
10.1±0.8 



kT (keV) r X 2 /dof 

1.57±0.01 ... 12666/189 

4.79±0.60 ... 8789/189 

3.31±g;J| 8104/188 

3.00±0.14 231/180 

1.56±0.01 ... 6677/182 

4.46±0.10 ... 3433/182 

2.99±0.12 2810/181 

1.52±0.05 2.74±0.17 127/179 



Table 2 

Best-fitting model to the non-dip spectrum 



N h 


kT 


r E CO 


E Hne 


8.6±3.1 


1.32±0.17 


2.29±?;f^ -12 


a 4 +0.19 



Table 3 
Best fits to the six spectra 



Spectrum N§ B N% PL N° PL f X 2 /dof 





Non-dip 


8.6±3.1 


8.6±3.1 


0.0 


0.0 


54/45 


280- 


-300 c s" 1 


8.6±3.1 


8-6^ 




0.123 ±0.009 


19/45 


260- 


-280 c s" 1 


8.6±3.1 


8.6±I;t 


148±10 


0.300 ±0.012 


63/45 


220- 


-240 c s" 1 


oo + 15 
ZJ -12 


8.6±H 


120±6 


0.504 ±0.005 


53/44 


100- 


-120 c s' 1 


> 10 3 


51±2 


53±3 


0.793 ± 0.006 


44/45 


80- 


-100 c s" 1 


> 10 3 


154±5 


60±6 


0.828 ±0.006 


28/45 
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